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Abstract 

In massive gravity the so-far-found black hole solutions on Minkowski space 
happen to convert horizons into a certain type of singularities. Here we explore 
whether these singularities can be avoided if space-time is not asymptotically 
Minkowskian. We find an exact analytic black hole (BH) solution which evades 
the above problem by a transition at large scales to self-induced de Sitter (dS) 
space-time, with the curvature scale set by the graviton mass. This solution is 
similar to the ones discovered by Koyama, Niz and Tasinato, and by Nieuwen- 
huizen, but differs in detail. The solution demonstrates that in massive GR, 
in the Schwarzschild coordinate system, a BH metric has to be accompanied 
by the Stiickelberg fields with nontrivial backgrounds to prevent the horizons 
to convert into the singularities. We also find an analogous solution for a 
Reissner-Nordstrom BH on dS space. A limitation of our approach, is that 
we find the solutions only for specific values of the two free parameters of the 
theory, for which both the vector and scalar fluctuations loose their kinetic 
terms, however, we hope our solutions represent a broader class with better 
behaved perturbations. 



1 Introduction and Summary 



According to the representation theory of the Poincare group in 4D, a massive spin-2 
state has to have five degrees of freedom; these can be thought of as the helicity-O, 
±1, ±2 states. A good Lagrangian for the massive spin-2 has to be able to describe 
these states. The Fierz-Pauli mass term pQ is the only ghost- and tachyon-free term 
at the quadratic order that describes the above 5 states [2|. However, in the zero mass 
limit it does not recover the linearized Einstein's gravity, since the helicity-O mode 
couples to the trace of the matter stress-tensor with strength equal to that of the 
helicity-2; this is called the vDVZ discontinuity [3J. If true, it would rule out massive 
gravity on the grounds of solar system observations. However, Vainshtein |4] argued 
that the troublesome longitudinal mode is suppressed at measurable distances by 
nonlinear effects, making the nonlinear theory compatible with current empirical 
data [5j. On the other hand, in a broad class of models, the same nonlinear terms 
that are responsible for the above-mentioned suppression, give rise to an instability 
known as the Boulware-Deser (BD) ghost [6j. This ghost appears as a 6th degree 
of freedom in the theory, and even though it is infinitely heavy on the Minkowski 
background, it becomes sufficiently light on locally nontrivial backgrounds, thus 
invalidating the theory [7J El El [10]. 

More recently, however, using the effective field theory formalism of [TJ, it has 
been shown in Ref. that there exists a two parameter family of nonlinear 

generalization of the linear Fierz-Pauli theory, that is free of the BD ghost order-by 
order and to all orders, at least in the decoupling limit. 

Most importantly, it was shown in Ref. [T2J that the absence of the BD ghost in 
the decoupling limit is such a powerful requirement that it leads to the resummation 
of the entire infinite series of the terms in the effective Lagrangian. As a result, a 
candidate theory of massive General Relativity free of BD ghost, was proposed [12J. 

Using the Hamiltonian analysis in the unitary gauge it was shown that for a 
certain choice of the free parameters of the theory, and in the 4th order in nonlin- 
earities, the Hamiltonian constraint that forbids the BD ghost is maintained in the 
theory of [12] . Note that the quartic order is special, since the lapse necessarily en- 
ters nonlinearly in all massive theories precisely in this order [8] , and it may appear 
that the hamiltonian constraint should necessarily be lost then. In spite of this, the 
constraint is maintained in a subtle way for special theories, as was shown for a toy 
model in [11] . and shown in the 4th order for massive GR in |12] . 

The existence of the Hamiltonian constraint to all orders in the unitary gauge, 
and for generic values of the parameters, was shown in Ref. [13j, using the method 
of dealing with the lapse and shift proposed in |TT| [12"] . Moreover, Ref. [13] has 
also argued for a secondary constraint, that follows from the conservation of the 
Hamiltonian constraint 0- Very recently, the existence of the secondary constraint 
was explicitly confirmed in Ref. |15| . 

The absence of the BD ghost among the local fluctuations of the theory of [T2] 

1 The argument of the existence of the secondary constraint was challenged in |14) . 
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in a generic gauge has been shown using the Stiickelberg decomposition [16j, as well 
as the helicity decompositions [T7] to quartic orders in nonlinearities (in the latter 
two references, previous misconceptions in the literature claiming the presence of 
the BD ghost were also clarified). Motivated by the above developments, in the 
present work we will proceed to study certain subtle properties black holes (BH) in 
the theory of [12]. 

In the unitary gauge Lagrangian of the theory the object = g^ v — 77^, is 
the gravitational analog of the Proca field of massive electrodynamics, describing all 
the five modes of the graviton. The diffeomorphism invariance can be restored by 
introducing the four scalar fields a (the Stiickelberg fields) [T8l[7l[T9]. and replacing 
the Minkowski metric by the covariant tensor d tl (j) a 'd v 4> b r}ab 

9^ = d^ a d v <p b r lab + H^ (1) 

where H^ v denotes the covariantized metric perturbation, and rj a b = diag(— 1, 1, 1, 1). 
The existence of the 4 Stiickelberg scalars <p a in this theory leads to the existence 
of new invariants in addition to the ones usually encountered in GR (Ricci scalar, 
Ricci tensor square, Riemann tensor square, etc); one new basic invariant is I ab = 
g^d^dyfy . Note that the unitary gauge is set by the condition <p a = x^d®. In this 
gauge, I ab = g^d^S^. Hence, any inverse metric that has divergence (even those 
which are innocuous in GR) would exhibit a singularity in the invariant I ab . Is 
this singularity of any significance? The singularity in the above invariant does not 
necessarily affect the geodesic motion of external observers - the geodesic equation 
is identical to that of GR, and due to its covariance, one could remove from it what 
would have been a coordinate singularity in g^ u in GR. However, one would expect 
the singularities in I ab = g» v 5^5^ to be a problem for fluctuations around classical 
solutions exhibiting it. Since g^ u could change signs on either side of the singularity, 
this could lead to emergence of ghosts and/or tachyons in the fluctuations around 
a given classical solution. In what follows, we will take a conservative point of view 
and will only accept solutions that have non-singular I ab . These arguments, in a 
somewhat different form, have already been emphasized recently by Deffayet and 
Jacobson |2"U] . 

The above arguments give rise to the following seeming puzzle. On the one 
hand, according to the Vainshtein mechanism [4J , spherically symmetric solutions of 
massive gravity should approximate those of GR better and better, as we increase 
the mass of the source and come closer to it. This would imply that the metric of a 
BH near its horizon should very much be similar to that of GR. On the other hand, 
the conventional Schwarzschild metric - if it were the solution of massive gravity in 
unitary gauge - would be singular according to the arguments above. 

We reiterate this central point in more general terms: In order for a metric 
to qualify as a valid description of a BH configuration, the physical singularities 
must be absent at the horizon. Then, in the unitary gauge of massive gravity the 
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Schwarzschild-like metric 



ds 2 = -(1 - f)dt 2 + 



dr 2 



+ r 2 dQ 2 , with e.g. / = r g /r, 



(2) 



1-f 



cannot be a legitimate BH solution of the theory. The same applies to the metric of 
de Sitter (dS) space in the static coordinates for which / = m 2 r 2 . 

Recently, interesting BH solutions of massive gravity have been found in Refs. 
|21[ |22| [23] (for other interesting solutions, which will not be discussed here, see, 
|24] - |32]). Following Koyama, Niz and Tasinato (KNT) |21J . one can start in the 
unitary gauge, and consider a most general stationary spherically symmetric metric. 
Then, using the method developed by KNT, very interesting full non-linear solutions 
for stars and black holes with Minkowskian asymptotics were found by Gruzinov 
and Mirbabayi in [23]. These solutions do exhibit the Vainshtein mechanism, and 
therefore are potentially viable classical solutions for stars and other compact objects 
in massive gravity (although their stability still remains to be studied). Nevertheless, 
it is not clear, as emphasized in [23], whether these are appropriate solutions for BHs. 
Even in the best case solution, when all the GR invariants are finite, the invariant 
g^dn^dvrfPrjab diverges, [23]. As noted above, this divergence does not affect the 
geodesic motion of any external observer, however, we expect it to be a problem for 
fluctuations. 

Could there be any solution that avoids the above issue? The answer is positive, 
and the resolution is in the identification of the unitary gauge to the coordinate 
system in which the black hole has no horizon (the Kruskal-Szekeres, Eddington- 
Finkelstein, or Gullstrand-Painleve systems come to our mind). The most convenient 
one for our purposes is the Gullstrand-Painleve (GP) system, in which the metric 
has the following form 



and is free of horizon singularities. It corresponds to the frame of an in-falling 
observer and covers half the whole space (for either choice of sign). 

Furthermore, if one has the metric Q as a solution in unitary gauge, then 
the coordinate transformation to the metric ([2]) will lead us to a background with 
cj) a 7^ x a . This means that if the configuration is described by the metric (J2]), the 
presence of a halo of helicity ±1 and/or fields around the BH is unavoidable. We 
will show in the present work that massive gravity in unitary gauge admits BH 
solutions precisely of this type. 

Interestingly, the dS-Schwarzschild solution found in [2Tj [33] do happen to satisfy 
our conservative criterion of non-singularity. However, the solution that we present 
here is not among the ones of [21, 33j. 

One more point worth emphasizing is that the BH solutions of [23] do exhibit 
the "helicity-0 hair" (e.g., produce an extra scalar force), while the ones found in 
[2Tj [33] . and in the present work do not. The status of the "no-hair" theorems in 




(3) 
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GR with the galileon field (which should capture some properties of the helicity-0 
of massive gravity) will be discussed in Ref. [34J. 

A limitation of our work is that we only manage to find these exact analytic 
solutions for a specific choice of the two free parameters of massive gravity. Such 
a choice is peculiar since on the obtained background, as we will show, the kinetic 
terms for both the vector and scalar fluctuations vanish in the decoupling limit. 

This fact would imply infinitely strong interactions for these modes (unless these 
modes happen to be nondynamical to all orders, e.g., due to the specific choice of 
the coefficients of the theory). Because of this issue, we would like to regard the 
solutions obtained here as just examples demonstrating how non-singular solutions 
should emerge. We also hope that our solutions are representative of a broader class 
of solutions which may have better behaved fluctuations. 

In this regard, there seems to be a few directions in which the studies of massive 
gravity BH's can be extended. First, one could look at the metric in the unitary 
gauge which would be some generalization of the Kruskal-Szekeres form. Second, 
one can extend the massive theory of [12] by adding more degrees of freedom to 
the existing 5 helicity states of massive graviton. In fact, two consistent extensions 
have already been discussed so far: (I) adding one real scalar field that makes 
the graviton mass dynamical [28J; (II) adding one massless tensor field with two 
degrees of freedom [35j that makes the internal space metric of the Stuckelberg field 
dynamical (bigravity). In the latter case cosmological solutions were found recently 
in [31] and [32], while BH's were studied in [36]. 

The work is organized as follows. Section 2 gives a brief review of the theory 
of massive gravity [12J. In section 3 we find an exact Schwarzschild-de Sitter solu- 
tion, and in section 4 an exact Reissner-Nordstrom-de Sitter, solution which have 
nonsingular I ab . These solution are similar to those discovered by Koyama, Niz and 
Tasinato, and by Th. Nieuwenhuizen, but differ in detail: our dS solution has no 
ghost even though the vector field is present (compare to [33]). Moreover, on the 
obtained solution the singularities in the invariant I ab are absent (compare to [22J). 
In the Appendix A we give another exact Schwarzschild solution that asymptotes 
to a conformaly rescaled Minkowski space, and briefly mention its peculiarities. In 
the Appendix B we discuss fluctuations on the selfaccelerated solution of section 3. 



A massive graviton is described by the Lagrangian density of [12] specified below 



where U is the potential for the graviton that depends on two free parameters 0^4 



2 The Theory 




(4) 



U(g,<f> a ) = (U2 + a 3 U 3 + a A U±) 



(5) 
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and the individual terms in the potential are defined as follows: 



U 2 = [/C] 2 -[/C 2 ], (6) 
U 3 = [/C] 3 -3[/C][/C 2 ] + 2[/C 3 ], (7) 
U A = [/C] 4 -6[/C 2 ][/C] 2 + 8[/C 3 ][/C]+3[/C 2 ] 2 -6[/C 4 ], (8) 

where K.£(g, <p a ) = 5^ — a/ g m da^ a dv^ h ^ab'i rectangular brackets denote traces, [K] = 
Tr(/C) = K,^. The above potential is unique - no further polynomial terms can be 
added to the action without introducing the BD ghost. 

The tensor H^ u represents the covariantized metric perturbation, as discussed 
in the introduction, which reduces to the h^ v in unitary gauge. While in a gauge 
unfixed theory we have 

H^ u = 9^u - d^ a d u (j) b ri ab . (9) 

Moreover, U is constructed in such a way that the theory admits the Minkowski 
background 

g,u = t^, p = x»6«. (10) 

Hence, it is natural to split </>'s as the background plus the 'pion' contribution 
cf) a = x a — Tc a , and as it was already mentioned in the introduction, the unitary 
gauge is defined by the condition ir a = 0. In the non-unitary gauge, on the other 
hand, it proves to be useful to adopt the following decomposition 

mA a + fry 

11 = J3 ' C 11 ) 

where A^ 1 describes in the decoupling limit the helicity ±1, while 7r is the longitudinal 
mode of the graviton (in the decoupling limit [7], M v \ — > oo and m — > 0, while 
A 3 = M p \m 2 is held fixed). This limit captures the approximation in which the 
energy scale is much greater than the graviton mass scale, E ^> m. 

For convenience, in what follows, we define the coefficients a and f3 which are 
related to those of §5§ by 013 = —(—a + l)/3 and «4 = — /3/2 + (—a + 1)/12. 
For generic values of the parameters a and (3 the theory exhibits the Vainshtein 
mechanism, as show in the decoupling limit and beyond [2T| [26] . As was 
emphasized in [TT], for one special choice, a = (3 = 0, the nonlinear interactions 
vanish in the decoupling limit with fixed A, leaving the theory weakly coupled (i.e., 
no Vainshtein mechanism) in this limit. For this particular choice of the coefficients 
the action of massive gravity with the potential ([6])- ([8]) (which can be rewritten in 
terms of just [K] and tuned to it cosmological constant [25], referred as a minimal 
model in Ref. [25]), was shown not to exhibit the Vainshtein mechanism also away 
from the decoupling limit [21] . 



5 



3 A Black Hole on de Sitter 



In this section we present the Schwarzschild-de Sitter solution in the theory of mas- 
sive gravity described above. The obtained solution is free of singularities (except 
from the conventional one appearing in GR). 

For convenience we choose unitary gauge for the metric. In this gauge the sym- 
metric tensor g^ u is an observable describing all the five degrees of freedom of a 
massive graviton. The equations of motion in empty space read as follows 

+ m 2 X^ u = , (12) 

where X^ v is the effective energy-momentum tensor due to the graviton mass, 

X(t» = ~\ £flV -JCw + at (k 2 ^ - KK,^ + ([/C] 2 - [/C 2 ])^j (13) 
+6/3 (lC% - KK\ V + \K„ V ([/C] 2 - [K 2 ]) - l 9lw ([/C] 3 - 3[/C][/C 2 ] + 2[/C 3 ])) 

Using the Bianchi identities, from ( TT2l) we obtain the following constraint on the 
metric 

m 2 V^ = 0, (14) 

where V M denotes the covariant derivative. 

In order to obtain the expression for we make use of the fact that the La- 
grangian is written as the trace of the polynomial of the matrix /C£. Thus, following 
the method by Koyama, Niz and Tasinato [21j, we choose the basis which diago- 
nalizes the expression appearing under the square root in the definition of /C^ [One 
should bear in mind that this is not a coordinate transformation, but rather a trick 
to simplify the procedure of getting the equations of motion]. As a result, the po- 
tential becomes a function of the components of the inverse metric, rather than the 
combination of square roots of matrices. Having done this, one is free to vary the 
action with respect to the inverse metric components to obtain explicit expression 
for (fT2|) . Since these expressions are quite cumbersome we will not give them here. 

Below, we concentrate on one particular family of the ghost-free theory of massive 
gravity in which there is the following relations between the two free coefficients: 

a 2 

/3 = -y (15) 

That this choice of the coefficients is special was first shown by Th. Nieuwenhuizen 
[22] (see also |23]). In particular, it was shown in [22J that for this choice the equation 
(JHJ) is automatically satisfied for a certain diagonal (in spherical coordinates) and 
time-independent metrics. It is interesting, however, that the above property persists 
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for a more general class of non-diagonal spherically-symmetric metrics written as 
follows: 

ds 2 = -A(r)dt 2 + 2B(r)dtdr + C{r)dr 2 + w 2 r 2 dQ 2 , (16) 

where w is a constant, while A(r), B(r) and C(r) are arbitrary functions. 

In subsection 3.1 we find an exact de Sitter solution to f|T2|) . and in subsection 
3.2 we find an exact BH solution on the obtained dS background. Note that the dS 
background is entirely due to the graviton mass. 

3.1 The de Sitter Solution 

We note that we would find an exact dS solution if we required that 

m 2 X^ = Xg^, (17) 

where A is some constant. The solution of the equations ( JT2|) that also satisfies ( |T7|) 
with a positive but otherwise arbitrary a is given by 

ds 2 = -n 2 dt 2 +| - dr ± k\ — , a r mrdt j + , - ^ r 2 dVL 2 . (18) 
ya + 1 V3a(a + 1) J (a + l) 2 K ' 

Here, k is a positive integration constant. It is straightforward to check that for ( fT8|) 
we have X^ v = (2/a) gr^, leading to the expression for the Ricci scalar 

R=-m 2 , (19) 

a 

as expected. Hence, this is a dS space with curvature scale set by the graviton 
mass and one free parameter a. One could imagine that m ~ (0.1 — 1)Hq, and 
a ~ (0.01 — 1), in which case the obtained dS solution (if stable) could describe 
dark energy. 

Up to a rescaling of the coordinates, the expression ( jl8p looks exactly like the de 
Sitter solution of GR written in the Gullstrand-Painleve frame. Either ± solution 
covers half of dS space. One can rotate the obtained solution to the static coordinate 
system at the expense of nonzero Stiickelberg fields. This will be done in the next 
subsection. In either form, the solution has no additional singularities. 

3.2 Schwarzschild-de Sitter Background 

Having the solution of the previous subsection worked out, it is straightforward to 
show that the system of equations f[T2"j) admits the following exact solution 

ds 2 = -K 2 dt 2 +\&dr± k\ ^ + —m 2 r 2 dt ) + a 2 r 2 dVi 2 , (20) 
\ V ar 3a 
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where a = a/(a + 1), and as before, k is an integration constant. In order to bring 
this solution to a more familiar form let us perform the following rescaling 



a + 1 
r — > r, 

a 

dt — > —dt. 



The resulting metric reads 



(21) 



ds 2 = -dt 2 +[dr±\l r -^ + —m 2 r 2 dt \ + r 2 dtt 2 . 

r 3a 



(22) 



This is the Schwarzschild-de Sitter solution in the GP coordinates. 
However, the above rescaling takes us away from the unitary gauge 

1 



t I 1 It 

K 

r — > r H r. 

a 



(23) 
(24) 



In terms of the 'pions', 71^ = x M — M , which can be decomposed as 7r M = (mA M + 
<9 M 7r)/A 3 , we have 



7T 



A 3 

Y 

0. 



(25) 
(26) 



The fields in f[26|) correspond to the canonically normalized fields carrying the helicity 
eigenstates in the decoupling limit. 

Let us now rewrite our solution into a more familiar coordinate system. The 
metric can be transformed to a static slicing by means of the following coordinate 
transformation 



dt^dt + f'(r)dr, 



(27) 



with f(r) = -gai/goo given by 



f'(r) 



3a 



1 — — — ^-m 2 r 2 

r ia 



(2* 



The resulting expression for the metric reads as follows: 
ds 2 = - 



l- r JL 
r 



* r ia 



+ r 2 dQ 2 



(29) 
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This is nothing but the metric of the Schwarzschild-de Sitter solution of GR in the 
static coordinates. However, this metric should be accompanied by a nontrivial 
backgrounds for the Stiickelberg fields. Indeed, it is evident that (12 7p gives rise to 
the shift 5(f) = f(r). In turn, this gives rise to a background for the 'vector mode' 

A = -—f(r), (30) 
Km 

A 1 = 0. (31) 

This particular field assignment has been chosen according to scaling in the decou- 
pling limit. Namely, f(r) vanishes in the decoupling limit linearly in m hence it was 
ascribed to the "vector mode". 

We would like to make two important comments in the remainder of this section. 
The first one concerns the integration constant k. Although, all the invariants of 
GR are independent of k, the new invariant that is characteristic of massive gravity 

I ab = g^d^ a dA\ (32) 

does depend on this integration constant; in the unitary gauge I ab is just the inverse 
of the GP metric (fT8|) which reads as follows 

/ -\ ± 1 o+l /o+lrs ~ 2 o \ 

ft K a y a r 3a (a+l) z 

± 1 a±l feUli + 2 a* ^ (£±l) 2 (l_2±li_ 2 a^ m 2 r 2\ Q 
k, a \J a r oa (a + l) z V a / V a r 3a (a-\-l) z J 

V o o mXxW 



Thus, the backgrounds with different values of k correspond to distinct super- 
selection sectors labeled by the values of I ab . 

The second comment concerns the issue of small fluctuations on top of this 
solution. One may worry that the scalar perturbations on this background are 
infinitely strongly coupled in the light of the results of |24j. In the latter work it 
was found that, for the parameters chosen as in 015p . the de Sitter background has 
infinitely strongly coupled fluctuations in the decoupling limit. However, we should 
point out that the self-accelerated background discussed in this section is different 
from that studied in [23]. This distinction is manifest in (13"T|) by the presence of the 
background for Aq, which vanishes in the case of [24J. 

Still, one could argue that it is unnecessary to perform the transformation of 
variables (|27|) responsible for this difference, and limit oneself to the rescaling of the 
coordinates 

a + 1 1 I \ 

r, t->-t, (33) 

a k 

which clearly does not give rise to the vector background. As a result, the 'pion' 
configuration will become similar to that of |24| . while the metric itself will be quite 
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different, namelj{§ 

ds 2 = -(1 - m 2 r 2 )dt 2 + 2mrdtdr + dr 2 + r 2 dVt 2 . (34) 

Now, if we were to take this metric as the one in which the decoupling limit should 
be taken, then we would find that the gauge freedom that is left in this limit 

9pv-* 9nv+ d (n£v)i ( 35 ) 

would not be enough for bringing (134"]) to the form of de Sitter space in either 
conformal or static slicing. Furthermore, the canonically normalized (13~4|) diverges 
in the decoupling limit, in such a way that this divergence can be isolated only in 
the vector mode. If so, then no conclusion can be drawn about the perturbations 
around our background based on the results of [23]. This, on the other hand, does 
not necessarily imply that the fluctuations are fine. As we show in the Appendix B 
the vector and scalar fluctuations may be infinitely strongly coupled. 



3.3 From Gullstrand-Painleve to Kruskal-Szekeres 

In this section we ask the question whether the BH solution of the GP form could 
be analytically continued to cover the other half of the space-time as well. This 
can be done by going to the Kruskal-Szekeres (KS) coordinates and analyzing the 
Stiickelberg fields. 

Let us start addressing this point by considering the following background 

ds 2 = -dt 2 GP + (dr + ^-dt G p^j +r 2 dQ 2 , <f) a = x a . (36) 

First we rewrite the metric in static slicing by performing the following change of 
the time variable 

t GP = t + 2r g J^- + r 3 ln ( V ^ | ■ (37) 



As a result the metric takes on the Schwarzschild form. In order to go to KS 
coordinates we use reparametrizations identical to the one used in GR 

— -l) e r/r ° =X 2 - T 2 , (38) 

r 9 J 

t = r 9 \n (|±T) . (39) 

For the analysis of the 0's in KS coordinates we make the 'near the horizon' ap- 
proximation r/r g — > 1, since this is the the region of our interest. In this limit 



2 For simplicity we set r„ = and drop the numerical factors. 
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the above coordinate transformation simplifies to (near the horizons T = ±X, with 
signs corresponding to the black- and white-hole respectively) 

— - 1] = -(X 2 -T 2 ), (40) 
r g J e 

t = r g ln . (41) 

As a result the 0's take the following form (using the fact that T 2 - X 2 is small) 

0° = 2r g lnpf + T) + r s (ln(l/4) + 1), (42) 
f = r g (l + - e (X 2 -T 2 )^. (43) 

Notice that 0° is singular on the horizon of the white hole (while being regular on 
the black hole horizon). The metric in these coordinates is given by 

ds 2 = \ r f-^-L{-dT 2 + dX 2 )+r 2 dQ 2 (44) 

The invariant I ab = g^d^d^ on the above background is singular at X = — T, 
corresponding to the horizon of the white hole. 

If one takes the original GP metric (I3"6"j) to describe the white hole instead of 
the black hole (this is achieved by flipping the relative sign of the expressions in 
parentheses) then after analytical continuation to KS coordinates the singularity 
will appear on the black hole horizon rather than on the one of the white hole. The 
generalization of this arguments for the case of the dS is straightforward. 



4 Reissner-Nordstrom solution on de Sitter 

The ghost-free theory of massive gravity with /3 = —a 2 /6, upon its coupling to the 
Maxwell's theory of electromagnetism, possesses the following Reissner-Nordstrom 
solution on dS space 

ds 2 = -dt 2 + I adr ± \l ^ + — m 2 r 2 - -^dt I + aVdft 2 , (45) 

ar 3a orr 1 1 

with a = a /(a + 1) and the electromagnetic field given by 

E=Q and B = 0. (46) 

In order to normalize the radial coordinate appropriately and to rewrite the solution 
in the static slicing, one needs to perform the rescaling 

r 

(47) 

a 
11 



supplemented with the following transformation of time 



ill- ->J-/--h/(>) ( //-. with /'(,) = -— = ± V r +3 " m2r2 ^ . ( 48 ) 



£00 1 - ^ - fmV 2 + ^ 



As a result the metric takes the familiar form 

<U = ( 1 - ^ - ^mV + ^)^ 2 + — — + r 2 dQ 2 , ( 19) 



r 3a 

while the Stiickelberg fields become 



0° = * + /(r), (50) 

4> T = r + -r. (51) 
a 

In this reference frame the electromagnetic field is 

E = -Q- and B = 0. (52) 
And, for obvious reasons the actual charge should be defined by Q = Q/a. 
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A Schwarzschild-like Solution 

In this appendix we concentrate on a different choice of the parameters 

a 2 

n = - T , (i) 

for which some exact solutions can also be obtained. In particular, we show that 
there exists an exact non-singular and asymptotically flat BH solution. We choose 
the unitary gauge and find that X tt and X tr from eq. f|T2|) vanish identically on the 
ansatz (fT6|) . for w = a /(a + 2). Suggesting, that there exist fluctuations which are 
infinitely strongly coupled. 
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Then, the solution to the full set of the equations of motion takes the following 
form 



(a + 2Ya 



3^,2 



± 



+ 



(a + 2) 5 + a 5 5 

2a{a + 2) 
a + 2f + a 5 5 

,2 

' 1 + 



1 - 



r g {a + 2) 



rot 



dt 2 



(II) 



cr 



a + 2) 



r g {a + 2f 



ar((a + 2) 5 + a 5 5) 



dr 2 



a 2 r 2 



^dn 2 , 

(a + 2) 2 



where 5 and r g are positive integration constants. The transformation to the Schwarzschild 
coordinates is carried out in a way similar to the previous section 



r — > 



a + 2 



a 



dt -> £(dt + f'(r)dr), with ( 2 

As a result the metric takes the conventional form 

dr 2 



(a + 2) 5 



+ 



a' 



[a + 2} 



ds' 



1 - 



dr + 



l - 



+ r W 



however, this should be accompanied by the 'pion' configuration 



7T 

A° 
A 1 



A3 



[l-C)t 2 --r 2 
a 



A 3 

-C-/W 
m 

0. 



(HI) 



(IV) 



(V) 

(VI) 
(VII) 



It is interesting that there exists a choice of parameters for which one gets the 
background metric identical to that of GR, supplemented with the 'pion' fields listed 
above. All the invariants are regular on this solution (away from the singularity in 
the center). It should be pointed out that ([Tl]) is the only background among those 
given in [21] with vanishing cosmological constant. This solution, however, has 
infinitely strongly coupled fluctuations; this and related issues will be discussed in 
a forthcoming paper |38| . 



B Fluctuations 

In this section we study the fluctuations on the backgrounds of section 3. The 
analysis is done in the decoupling limit [7] 

m ->• 0, M p i -> oo, with A 3 = M pl m 2 - fixed. (VIII) 
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In this limit the Lagrangian can be decomposed into tow pieces. The first describes 
the dynamics of the helicity-O, ±2 modes and their interactions with each other. 
While the second one accounts for the helicity-±l modes and their nonlinear cou- 
plings to the helicity-O degree of freedom. 

The scalar-tensor Lagrangian, with the condition /3 = —a 2 /6, is given by [11] 

£st = -\h^h a , + hr (xg) - ^ - ^xg) . (ix) 

Here, the first term represents the linearized Einstein-Hilbert Lagrangian, while X's 
are defined in terms of the longitudinal mode as follows 

2 ^papcr^uf5pa^-af5 1 
^pap^^ufiaS^-afl^-pa^-jSi 

with 11^2, = d^d v ix and all the repeated indices contructed by the flat space metric. 

The scalar-vector Lagrangian, on the other hand, contains an infinite number of 
terms and schematically is given by 

£^-^ + E^§)", (X) 

n=l ^ ' 

where F^ v denotes the field strength of the helicity-1 mode, A^, and the whole La- 
grangian is invariant under the U(l) gauge transformation 5A^ = d^a. Remarkably, 
this expression has been recently resummed for the spherically symmetric ansatz 
[33j, making the analysis of the vector fluctuations possible. 

After expanding Cst, and the resummed version of £$v (eq. (C.5) of Ref. [33] 
which is too lengthy to be reproduced here) to the second order in perturbations 
around the backgrounds of section 3, we find that the kinetic terms for both helicity- 
and helicity-±l fields vanish identically, on the dS as well as the Schwarzschild-dS 
space. This does imply that these modes are infinitely strongly coupled in this limit 
(unless of course they are rendered nondynamical to all orders by some symmetry 
or constraint). Whether this problem can be remedied by going beyond the decou- 
pling limit, and/or by invoking kinetic terms due to quantum loops (which will be 
generated as long as they're not prohibited by symmetries), remains to be seen. 
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